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Abstract 

We discuss the restructuring of the BPS spectrum which occurs on certain 
submanifolds of the moduli/parameter space - the curves of the marginal 
stability (CMS) - using quasiclassical methods. We argue that in general 
a 'composite' BPS soliton swells in coordinate space as one approaches the 
CMS and that, as a bound state of two 'primary' solitons, its dynamics in this 
region is determined by non-relativistic supersymmetric quantum mechanics. 
Near the CMS the bound state has a wave function which is highly spread 
out. Precisely on the CMS the bound state level reaches the continuum, 
the composite state delocalizes in coordinate space, and restructuring of the 
spectrum can occur. We present a detailed analysis of this behavior in a 
two-dimensional J\f= 2 Wess-Zumino model with two chiral fields, and then 
discuss how it arises in the context of 'composite' dyons near weak coupling 
CMS curves in J\f= 2 supersymmetric gauge theories. We also consider cases 
where some states become massless on the CMS. 
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I. INTRODUCTION 



Centrally extended supersymmetry algebras admit a special class of massive representa- 
tions which preserve some fraction of the supersymmetry of the vacuum, and consequently 
form multiplets which are smaller than a generic massive representation. The states lying in 
these shortened (or BPS, after Bogomol'nyi-Prasad-Sommerfield) multiplets are extremely 
useful probes of the theory because on one hand their spectrum is determined almost entirely 
by kinematical constraints (i.e. by the central charges), while on the other the multiplet 
structure ensures their generic stability. More precisely, the fact that BPS states lie in short 
multiplets implies that they must remain BPS, unless a degeneracy of several BPS multiplets 
is achieved which can then combine to form a generic massive multiplet. In the absence of 
such an exotic scenario, the dynamics of the BPS sector forms a closed subsystem. 

The stability of BPS (particle) states follows from the fact that their masses are deter- 
mined by the superalgebra to be the expectation values of the central charge (s) Z^. Since 
the central charges are additive, this implies via the triangle inequality that a BPS state 
whose mass is 



Even at points where this equality is saturated there is no phase space for a physical decay. 
Thus one concludes that BPS particles are indeed stable. 

However, restructuring of the spectrum is nonetheless possible because of the existence 
of special submanifolds of the moduli/parameter space where the inequality (2) is saturated. 
Specifically, this allows for discontinuities of the spectrum with respect to changes in these 
moduli. Such changes are 'unphysical' in the sense that one shifts between different supers- 
election sectors. Nonetheless, one is often interested in considering such an evolution, as it 
may correspond to the extrapolation from a weakly coupled to a strongly coupled regime. 
In this case, the stability of BPS states can often be used to infer information about the 
strongly coupled region. The caveat of course is that one should not cross a submanifold 
where the bound (2) is saturated, and where restructuring of the spectrum may occur and 
BPS states may for example disappear. Such submanifolds are consequently known as curves 
of marginal stability (CMS), although their actual co-dimension in the moduli/parameter 
space will vary. 

Marginal stability curves, and the corresponding discontinuities of the BPS spectrum, 
are quite ubiquitous in theories with centrally extended supersymmetry algebras. Examples 
include: the existence of a CMS for the BPS soliton spectrum in general classes of two 
dimensional models discussed by Cecotti and Vafa [1]; and the CMS for the BPS particle 
spectrum in M = 2 supersymmetric gauge theories [2-5]. In the latter case an explicit 
demonstration of the discontinuity of the spectrum across these curves in the vacuum moduli 
space was provided in generic SU(2) theories by Bilal and Ferrari [4,5]. The realization of 
these dyonic states in terms of Type IIB string junctions has also led to the appearance of 




(1) 



is stable with respect to decay into BPS 'constituents' with masses Mi = \(Zi)\, 



(2) 
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marginal stability conditions in this context [6,7]. Furthermore, a discontinuity in the BPS 
spectrum of wall solutions in a Wess-Zumino model with the Taylor- Veneziano-Yankielowicz 
superpotential [8], which leads to a glued potential, was also observed recently by Smilga 
and Veselov [9,10]. The discontinuity arises in this function of the mass parameter 

- a feature also observed in some of the models to be considered in this paper. Finally, 
we also mention that marginal stability conditions have more recently been studied in the 
context of string compactification on manifolds with nontrivial cycles [11]. 

In order to illustrate the discussion with a particular example, we recall that the notion 
of marginal stability arises, in particular, in the Seiberg-Witten solution [2] of M— 2 su- 
persymmetric (SUSY) gauge theories (see e.g. [4]). In the simplest example of pure super 
Yang-Mills with gauge group SU(2), there is a one dimensional elliptical curve of marginal 
stability in the moduli space (see Fig. I). 




FIG. 1. A schematic representation of the moduli space for J\f= 2 SYM with gauge group SU(2) in 
terms of the vev u = (tr (f> 2 ) of the adjoint scalar 0. The W bosons only exist outside the shaded region, 
which consequently determines their stability domain. 

On crossing this curve by varying the moduli a restructuring of the spectrum of BPS 
states takes place. For instance, the electrically charged vector bosons only exist outside 
the CMS, and disappear from the spectrum in the interior region. To make these notions 
a little more general, we can define a "stability domain" as a submanifold of the moduli 
space in which a particular BPS state exists. This domain will always be bounded by a 
CMS. In this example, the W boson has a stability domain in the exterior region illustrated 
in Fig. I. On crossing the CMS from the stability domain, it is usually stated that the W 
bosons "decay" into a two particle state consisting of a monopole and a dyon with unit 
electric charge. This interpretation is a little awkward because for a particle to properly 
decay it must exist in the spectrum, at least as a quasi- stationary state, and this is not true 
after crossing the CMS. The question then arises as to exactly what happens on the CMS 
resulting in the apparent discontinuity of the BPS spectrum. 

In this paper we will suggest a physical interpretation of this phenomenon, which we 
summarize below. For this purpose, its convenient to continue with the W boson example 
to make the ideas more concrete. However, one should bear in mind that this system is 
not directly accessible to the semi-classical techniques that are used in this paper, because 
the CMS curve lies at strong coupling. Nonetheless, we will argue that there are several 
constraints ensuring, at least qualitatively, the generality of this behavior. 



3 



Specifically, the emerging picture is that when the moduli approach the CMS, the W ± 
states swell in coordinate space as they become more weakly bound. Near the CMS, but still 
within the stability domain, one can interpret the W ± as a composite particle built from 
two primary constituents (a monopole and a dyon of electric charge one), whose interaction 
can be described by a nonrelativistic (super)potential depending on the relative separation, 
within the framework of supersymmetric quantum mechanics (SQM) [12]. As one approaches 
the CMS, the separation of the two primary constituents diverges, while the bound state 
level reaches the continuum (i.e. the binding energy vanishes). Further motion after crossing 
the CMS leads to an SQM superpotential which fails to exhibit a supersymmetric ground 
state separated by a gap from the continuum. Consequently, the ground state in the sector 
with unit electric charge is no longer the one-particle W boson state but rather a set of 
non-BPS two-particle states forming a 'long' multiplet. 

We will argue that this picture is the general situation for CMS curves associated with 
BPS particle states. Namely, whenever a discontinuity occurs in the BPS spectrum at a point 
in the parameter space, then certain BPS states delocalize in coordinate space. Indeed, the 
phenomenon of marginal stability of BPS states involves, by definition, the alignment of 
central charges of primary states in such a way that the binding energy vanishes. In this 
context it is quite natural that crossing the CMS involves infrared effects, and the 'size' of 
the marginally stable state should diverge as the CMS is approached. However, within this 
general picture of derealization one can identify several different mechanisms underlying 
this behavior. 

The features are somewhat dimension dependent, so its convenient to focus first on 1+1D 
which will be our primary concern in this paper. We will then remark on certain aspects 
which distinguish the behavior in 3+1D in particular. Moreover, for solitons in 1+1D its 
convenient to distinguish two derealization mechanisms. 

(1) The first is when there are no massless fields relevant to the problem, and consequently 
one can describe the interactions of the primary constituents using non-relativistic collective 
coordinate dynamics with linearly realized supersymmetry and short range potentials. For 
a large class of systems (including the ones to be considered here), its possible to limit 
attention to just one collective coordinate - the relative separation r of the primary solitons. 
We then observe two characteristic dynamical scenarios for the near CMS dynamics: 

• In the first case, the short range potential is of deuteron-type which remains finite on 
the CMS but possesses a single bound state, whose wavefunction spreads out as the 
CMS is approached, while the level reaches the continuum at this point. On crossing 
the CMS, there is no longer a supersymmetric ground state reflecting the fact that the 
BPS state has disappeared from the spectrum to be replaced as a ground state by a 
non-BPS two-particle state with the same quantum numbers. 

• In the second scenario, the relative separation becomes an exact modulus on the CMS, 
and the potential therefore vanishes at this point. In this case, the composite state still 
delocalizes as the CMS is approached as the wavefunction becomes highly spread out. 
The state is however highly quantum mechanical and has no classical analogue. In this 
case, we also observe that the potential may support (in general different) composite 
states on each side of the CMS. 
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We will study a two-field model exhibiting both these dynamical regimes in subsequent 
sections. 

(2) The second derealization mechanism arises when there are massless fields involved, 
these being either the primary states themselves, or the fields via which their interactions 
are mediated. 

• In situations where massive primary states interact via massless exchange, we shall 
argue in Sec. 2 that attractive Coulomb-like interactions between the primary con- 
stituents must vanish on the CMS as a consequence of the structure of the BPS mass 
spectrum. 

• The second mechanism arises when one or more of the primary states is massless. This 
scenario may be taken as a special case of (1) in that massless points arise generically 
as co-dimension one submanifolds on the CMS curve. In this case it is not possible 
to reduce the effective dynamics to non-relativistic quantum mechanics, and one must 
consider the effective theory of the massless state. We note that more exotic examples 
of this behavior may arise (in higher dimensions) at Argyres-Douglas points [13] in 
J\f = 2 SYM, or more generally at second order critical points in supersymmetric 
theories. 

Although we have framed this discussion in the context of 1+1D field theories many 
of the features apply also in higher dimensions. In particular, restructuring of the BPS 
spectrum via derealization is apparently a generic phenomenon. However, an important 
distinction between 1+1D and, say, 3+1D is that in 1+1D an arbitrarily small attraction is 
sufficient to form a bound state while in 3+1D this is not the case. For this reason long range 
forces play a special role in 3+1D, where Coulomb-like attractive potentials are required to 
form bound states at an arbitrarily small effective coupling. We will discuss this scenario 
in the form of composite dyons in J\f— 2 SYM. The general arguments outlined above will 
be presented in Sec. II, while particular examples of different scenarios will be discussed in 
subsequent sections. 

In this paper we focus first on class (1) and present an exhaustive study of a particular 
two dimensional Wess Zumino model [14] with M= 2 supersymmetry of the type considered 
previously [15] in a related context. This is a simple model which exhibits composite solitons 
(kinks) and a corresponding CMS curve accessible to quasiclassical techniques. Thus it 
serves as an ideal arena to study in detail the effective SQM which determines the presence 
or otherwise of the composite soliton. The model involves two weakly interacting chiral 
fields. In the decoupling limit there are 'primary' BPS kinks for each field, which when 
quantized lead to short BPS multiplets containing one bosonic and one fermionic state (plus 
antiparticles). There are also 'composite' solitons which are combinations of the primary 
configurations. 

Switching on an interaction between the fields we see that the primary BPS solitons 
exist throughout the parameter space, while the composite solitons exhibit a finite stability 
domain bounded by the CMS. We analyze in detail the effective SQM which exhibits the 
composite soliton as a supersymmetric bound state, and verify the behavior described earlier 
with regard to the approach to the CMS. Its worth noting that in this model the structure 
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of the stability domain is quite complex. In particular, there are different dynamical regimes 
depending on which part of the CMS curve is crossed. In most composite state only 

exists on one side of the CMS. However, there are regions in the parameter space where 
stability domains for particular composite states meet on a CMS, and consequently (in 
general different) composite states can exist on either side. In this latter regime we observe 
that the relative separation of the primary states becomes a modulus on the CMS, as the 
potential vanishes. Therefore this model exhibits both scenarios outlined in (1) above. 

The advantage of dealing with the Wess-Zumino model is that all the features of the 
non-relativistic quantum dynamics can be calculated analytically. In the vicinity of the 
CMS we obtain the explicit form for the (super)potential describing the interaction of the 
primary solitons and are able to track the form of the bound state wavefunction right onto 
the CMS. Moreover, certain qualitative aspects are apparently rather model independent 
due to the constraints imposed by the BPS spectrum. To investigate the situation in 3+1D 
we consider explicitly Af= 2 SYM with gauge group SU(3) which contains a spectrum of 
primary and composite monopole (dyon) solutions. The two 'primary' monopole solutions 
are embedded along each of the simple roots of the algebra. An embedding along the 
additional positive root leads to a 'composite' dyon which becomes marginally stable on 
a CMS which is present in the semiclassical region. The major difference between the 
monopole case in four dimensions in comparison to the two-field model in two dimensions is 
the presence of massless exchanges resulting in a long range Coulomb-like interaction, which 
can lead to bound states as noted above. Recently there has been considerable interest in 
this system [16-22], in part because the composite dyon configuration is an example of a 1/4- 
BPS state in A/ = 4 SYM. This work, which has centered on the moduli space formulation of 
the low energy dynamics, has resulted in the detailed form of the long range interaction. We 
observe that, in accord with the general expectations of Sec. II, the attractive component of 
the long range force (the term oc 1/r in the effective potential) vanishes on the CMS, while 
a repulsive component (oc 1/r 2 ) remains. There is no attraction on the other side of the 
CMS, the term oc 1/r change its sign. Thus, a BPS bound state which exists in the stability 
domain on one side of the CMS becomes more and more delocalized when approaching the 
CMS, and ceases to exist on the other side. Accounting for the fact that long range forces 
are now crucial, we observe that the qualitative picture is nonetheless quite similar to the 
two-field Wess-Zumino model, in that the composite state delocalizes on approach to the 
CMS. 

The layout of the paper is as follows: In Sec. II we present some general arguments 
constraining the dynamics of primary solitons near a CMS. Using these results we discuss, 
in a simplified setting, the underlying mechanism involved in restructuring the spectrum, in- 
troducing the necessary notation and definitions in passing. In this section, we also consider 
the embedding of the effective SQM superalgebra within the superalgebra of the field theory. 
As a specific example, we consider the realization of the M— 2 superalgebra with central 
charges in two dimensions in the two soliton sector. In this regard its worth remarking that 
this embedding shows explicitly how the presence of field theoretic central charges is crucial 
in allowing a linear realization of supersymmetry in the effective non-relativistic dynamics. 

Sec. Ill presents a detailed analysis of the (quasiclassical) solitons in the two-field Wess 
Zumino model with regard to their BPS properties. We calculate the form of the CMS, and 
prove that outside the stability domain the BPS solution corresponding to the composite 
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state ceases to exist. In Sec. IV we derive and discuss the SQM which describes the inter- 
action between the primary solitons in the vicinity of CMS and determines whether or not 
a supersymmetric bound state exists. We obtain analytic solutions for the superpotential 
and the bound state wavefunction. 

In Sec. V we consider the more complex situation of monopoles and dyons in A/ = 2 
SYM with gauge group SU(3), and review the form of the long range potential near the 
CMS [16-22]. The attractive component vanishes on the CMS, in agreement with the 
general arguments of Sec. II, while a repulsive component remains leading to derealization 
on the CMS even at the classical level. 

In Sec. VI we turn to the final derealization mechanism discussed in (2) above, which 
involves derealization due to a field becoming massless on the CMS. We consider a restric- 
tion of the two-field model, discussed in Sec. Ill which, when perturbed by a term which 
breaks H— 2 to Af= 1 supersymmetry, provides a simple example of this phenomenon. 

We collect some concluding remarks in Sec. VII, and discuss in particular the applicability 
of our results to marginal stability of the W boson, and also subtleties associated with 
extended BPS objects. 

II. SOLITON DYNAMICS NEAR THE CMS 

Before considering a specific model in detail, we first discuss some simple but quite general 
constraints which are useful in providing a qualitative guide to the dynamics appropriate to 
the near-CMS regime. 

A. Dominant interactions 

Consider the dynamics of two primary BPS solitons with masses M\ and M 2 near a 
CMS curve for the composite BPS soliton with mass Mi +2 . From the CMS condition 
that the binding energy vanishes, M 1+2 = Mi + M 2 , its clear that by going sufficiently 
close to the CMS, the relevant energy scales - kinetic and binding energy - can be made 
much smaller than the soliton masses. The system is then non-relativistic, and the effective 
dynamics is supersymmetric quantum mechanics on the space of collective coordinates of 
the configuration. With spherically symmetric interactions, the relevant part of this system 
can be reduced to one dimensional SQM associated with the relative separation r between 
the primary solitons. 

We can also deduce some generic features of the potential, in part from knowledge of the 
BPS mass spectrum. Firstly, its inconsistent for the potential to be of attractive Coulomb- 
like form on the CMS itself. This result follows straightforwardly from the incompatibility of 
the BPS mass spectrum with the structure of the bound state energy levels associated with 
a Coulomb-like potential. Indeed the quantum mechanical spectrum associated with the 
attractive 1/r potential will exhibit towers of closely spaced bound states, only the lowest 
of which can be BPS saturated. The Coulomb wavefunctions ip ~ r n e -r/n j ea j ^ Q ^q^^ 
state levels of the form e n oc —1/n 2 . In contrast, we know from the form of the BPS mass 
spectrum that on the CMS the lowest level in the tower must reach the continuum. Clearly 
the only way this can happen is if the 1/r attractive interaction vanishes on the CMS. 
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In other words, if attractive Coulomb-like forces are generically present, there must be a 
coefficient which we may identify as the distance to the CMS, 

V(r) const - (g* - /)_L + . . . , (3) 

where q is used to denote the appropriate charge and / is a certain function of the moduli 
equal to q 2 on the CMS. The ellipsis denotes higher order terms in 1/r. 

A second constraint is the requirement that the potential admits a normalizable bound 
state arbitrarily close to the CMS (inside the stability domain). This constraint is dimension 
dependent. While in 1+1D and 2+1D an arbitrarily small attraction can result in such 
bound state, this is not the case in higher dimensions. In 3+1D, in order to form a bound 
state the attraction must be strong enough, / drr (— V(r)) > h 2 /M. In particular, for the 
3+1D dynamics of dyons in SU(3) SYM, as we will see in Sec. V, the bound state is due 
to a Coulomb-like attraction at large distances [16-22]. Although according to Eq. (3) the 
effective Coulomb coupling diminishes on approach to the CMS, the bound state does exist 
even for an arbitrarily small coupling. 

In conclusion, from the simple arguments above we deduce that close to the CMS the 
dynamics is nonrelativistic and the long range component of the potential controlling the 
restructuring of the spectrum satisfies the following constraints. Firstly, on the CMS it either 
vanishes, or is repulsive. Secondly, the simplest way to form bound states in dimensions 
higher than 2+1D is for the potential to have an attractive long range form off the CMS. 

B. The restructuring mechanism 

To understand what happens to the spectrum in the near-CMS regime it will be useful to 
present a simple model which exhibits the relevant features. Specifically, we consider below 
the mechanism via which a restructuring of the spectrum can occur. 

Assume that the model under consideration contains a set of parameters (to be denoted 
generically as {//}), and admits BPS solitons at a certain value {/xo}- The parameters {/i} 
can be moduli, or some parameters in the action. The question is how can BPS solitons 
disappear from the spectrum under continuous variations of {/u}? Generally speaking, we 
would expect that if the BPS state exists at {ii }, it remains in the spectrum at least in 
some finite domain in the vicinity of {/U }. 

The argument is based on the multiplicity of the corresponding supermultiplet. Indeed, 
in the models to be considered below, the number of states in the BPS multiplet is twice 
smaller than the number of states in the non-BPS multiplet (this type of 'shortening' is 
typical). This means that if a BPS state is to become non-BPS, a factor of two jump in the 
number of states must occur. Generally speaking, this will not happen under continuous 
deformations of {//}, unless from the very beginning we had two BPS multiplets which 
become degenerate at a certain point in the parameter space and combine together to leave 
the BPS spectrum joint non-BPS multiplet. 

We are more interested in another scenario - when a BPS state becomes non-BPS at a 
certain critical point {//*}, without the pre-arranged doubling of the type mentioned above. 
Are we aware of any simple analogs of this phenomenon? 



8 



The answer is yes, a simple example has been known for a long time. We will discuss 
it here for two reasons: firstly, it nicely illustrates the generalities of the dynamical phe- 
nomenon discussed in the preceding subsection; and secondly, we will need to introduce the 
corresponding notation later anyway. The example can be found in supersymmetric quan- 
tum mechanics (SQM) with two supercharges introduced by Witten [12]. Consider a system 
(as motivated above) described by the Hamiltonian 



H 



1 



p 2 + (W') 2 + a 3 W"] , (4) 



where p = — id/dx, and W is a function of x with the prime denoting differentiation by 
x. Moreover, a% is the third Pauli matrix corresponding to the fact that [o"i,<r 2 ] forms an 
appropriate representation of the Grassmann bilinear. The function W will be referred to 
as the SQM superpotential. Two conserved (real) supercharges are 

Ql = 7! (p(Tl + w ' a2) ' Q2 = 72 {pa2 ~ w ' ai) ' (5) 

They form the following superalgebra, 

(Qi) 2 = (Q 2 ) 2 = # , {Qi,Q 2 } = 0. (6) 

If W has an odd number of zeros then the ground state of the system (4) is supersymmetric 
(i.e. the supercharges annihilate it) and unique. This is the analog of the BPS soliton. If 
W has no zeros or even number of zeros, the ground state is doubly degenerate and is not 
annihilated by the supercharges. The ground states in this case are analogs of non-BPS 
solitons. The unique versus doubly degenerate ground state in the problem (4) imitates 
"multiplet shortening". The transition from the first case to the second under continuous 
deformations of parameters is easy to visualize. Indeed, let us assume, for definiteness, that 

W(x) = In cosh a; — fix , W' = tanh x — \i . (7) 

At fj, = the derivative of the superpotential vanishes at the origin. As \x grows (remaining 
positive), the point where W vanishes shifts to the right, towards large positive values of x. 
The ground state wave function is supersymmetric and unique, 

fix 

As one approaches fi* — 1 from below this wave function becomes flatter on the right semi- 
axis; representing a swelling of the bound state in coordinate space. The corresponding 
scalar potential 

V(x) = - \(W'f - W" 
at pL — 0.98 and the ground state wave function are depicted in Fig. 2. 
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V|/(x) 



X 



V V(x) 

FIG. 2. The potential V(x) in the problem (4), (7) (solid line) and the corresponding ground state 
wave function (dashed line). The parameter fj, = 0.98. The units on the vertical and horizontal axes are 
arbitrary. 

The point // = 1 is critical. At n > 1 the wave function (8) at E — becomes non- 
normalizable, and the true ground state, coinciding with the continuum threshold, spectrum 
is doubly degenerate. The transition from one regime to another occurs through dereal- 
ization in that the zero of W, the equilibrium point xo, escapes to infinity. Note that 
dynamically the SQM problem under consideration is similar to that of deuterium. The 
potential well in Fig. 2 is at x < 1, but the tail of the wave function stretches very far to 
the right due to the fact that the E = level is very close to the continuum spectrum. 

We can make this somewhat more precise by introducing a "classicality parameter" £ 
defined as 



where xq is the classical minimum of the potential: W'(xq) = 0. The parameter £ may be 
interpreted as measuring the quantum correction to the curvature of the potential at the 
classical equilibrium point, i.e. the system is essentially classical if £ <C 1, while it is highly 
quantum if £ ^> 1. 

In the current example, we find that as we approach the critical point, 



and so the system indeed becomes highly quantum in this regime. In fact this feature is 
quite generic for short range potentials and may be viewed as an artifact of the remnant 
super symmetry. Specifically, since the mass term V F ~ W"(x) is linear in the superpotential, 
while the bosonic potential is quadratic, Vb ~ (W'(x)) 2 , for short range interactions the 
fermionic W" term in the superpotential (II B) will dominate for large separations. This is 
despite the fact that the fermionic term is a quantum effect (in field theory it corresponds 
to the 1-loop correction to the effective potential through integrating out the fermions). 
Thus, although the system becomes more and more weakly bound, in the CMS region the 
system enters a highly quantum regime where the classical minimum of the bosonic potential 
need not be relevant. Below we will see that exactly the same phenomenon occurs for BPS 
solitons near the CMS in a 1+1D Wess-Zumino model. 




(9) 
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C. Embedding of SQM within the field theory superalgebra 

To establish a link between the field-theoretical description of solitons on the one hand 
and the supersymmetric quantum mechanics of two nonrelativistic primary states on the 
other, we now consider the manner in which the quantum mechanical supercharges emerge 
from the full field-theoretical superalgebra. The fact that supersymmetry is realized linearly 
in the two soliton sector may be reinterpreted as the existence of a straightforward embed- 
ding of the SQM supercharges. Moreover, near the CMS the system becomes essentially 
nonrelativistic and we need keep only the leading term in an expansion in velocities. 

Although the arguments apply more generally, we consider for definiteness the realization 
of Af— 2 supersymmetry in two dimensions in the two soliton sector. Recall that the algebra 
contains four supercharges Q a , Q^ a (a = 1,2) and has the form [1,23] 

{Q a ,Ql} = 2(^°) a pP,, {Q a ,Qp} = 2i( 1 5 1 °UZ, {QlQ\} = 2i{ 1 ^) a pZ, (11) 

where P^ = (Pq, Pi) are the energy-momentum operators and Z is a complex central charge. 
We use the Majorana basis for 2x2 7-matrices, 

7 = a 2 , 7 1 = ^3, 7 5 = 7V = ■ (12) 

Modulo addition of the central charge, the algebra (11) can be viewed as a dimensional 
reduction of the M= 1 algebra in four dimensions. The SO(3,l) Lorentz symmetry in 3+1 
dimensions reduces in 1+1 to the product SO(l,l) xU(1)r where SO(l,l) is the Lorentz 
boost in 1+1 and U(l)# is a global symmetry associated with the fermion charge. More 
precisely, the Lorentz boost with parameter (3 acts on the supercharges Q a as follows 

Qp , (13) 

a/3 

while the U(1)a transformation with parameter 77 is 

Qp ■ (14) 

a/3 

Notice that the U(l)# transformations can be viewed as a complexification of the Lorentz 
boost (13). 

Its now convenient to introduce the Majorana supercharges Q l a , (i — 1,2; (Q l a Y = Q l a ) 
via the relation 

Q* = t ^{Q 1 a + i<&) , (15) 

where the phase factor Q~ ta / 2 contains an arbitrary parameter a, which we will fix momen- 
tarily. In terms of Q l a the algebra (11) has the form: 

{Qi, Q j p} = 2 W7 V P, + 2 i (7 5 7°W Z l > , (16) 

where the 2x2 real matrix of central charges Z l i is symmetric and traceless. It is related 
to the original complex Z as follows, 



Qc 



exp l-(3r 



Qc 



exp \ -rn 
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Ze~ ia = Z n -iZ 12 . (17) 
To consider representations of the algebra we use a Lorentz boost in 1+1 to put the 



system in the rest frame where Pi — > and P — ► M = yP M P^. Moreover, we can always 

choose the basis in XJ(1)r to put Z tj in the form Z tj = \Z\ . This amounts to fixing the 
phase a to be equal to the phase of the central charge, Z = \Z\ e ta . Then the algebra (16) 
takes the following component form, 

{Qlf = (Qif = M+\Z\, (Qlf = (Ql) 2 = M-\Z\, (18) 

with all other anticommutators vanishing, so that the algebra splits into two independent 
subalgebras. 

/,From (18) we see that \Z\ is a lower bound for the mass, M > \Z\. When M > \Z\ the 
irreducible representation has dimension four - two bosonic and two fermionic states. The 
BPS states saturate the lower bound, M B ps = \Z\, and in this case the second subalgebra 
becomes trivial and the representation is two-dimensional - one bosonic and one fermionic 
state [23]. 

How do all of these generalities help us with the problem of constructing the SQM near 
the CMS? In the vicinity of the CMS the difference M — \Z\ is small as compared to \Z\ 
and can be identified with the nonrelativistic Hamiltonian, 1 

H SQM = M-\Z\. (19) 

The second subalgebra in Eq. (18) with supercharges Q\ and Q\ then coincides with that 
of the standard SQM, see Eq. (6). In the first subalgebra the operator M + \Z\ can be 
substituted by 2\Z\ up to relativistic corrections. Consequently, the first subalgebra just 
leads to a generic multiplet structure (in this case just duplication) for every state found in 
the SQM. 

In Sec. IV we will find all the supercharges in the 1+1 example as explicit functions of the 
moduli from field-theoretic solutions for two solitons, u and v. Near the CMS, where their 
relative motion is nonrelativistic, the result can be compared with the quantum mechanical 
realization of the superalgebra (18). For -Psqm — M — \Z\ we take the expression which 
generalizes Eq. (4) to include a mass parameter, 



Hsqm — 



1 



2M, 



p 2 + (W) 2 + a 3 W"} , (20) 



where the superpotential W depends on the separation s = z u — z v , the conjugate momentum 
p = — id/ds, and M r is the reduced mass, 

M r = — — . 21 

M u + M v K } 



Note that we view M as a Hilbert space operator. 
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Then a realization of the superalgebra can be chosen in the form (cr* and 7* are two sets of 
Pauli matrices): 



q\ = J2\z\n ®a 3 , = / 2 | Z | T20a3 ^ 



Q2 



[p^ + W^s)^] , Qf = I 



y/2~W r 



\p<j 2 -W'(s) a,} 



(22) 



The realization (22) explicitly indicates a factorization of both the bosonic and fermionic 
degrees of freedom associated with the center of mass of the system. We can also include 
dependence on the total spatial momentum P 1 through a Lorentz boost (13) with tanh/5 = 



Px/y'M 2 + P 2 . 

A couple of comments are now in order. Firstly, its clear from this construction that 
the SQM can only be realized linearly in BPS sectors with a non-vanishing central charge. 
Otherwise, one has Q = y/~M ip (with vp a fermionic operator) in the nonrelativistic limit, 
implying a nonlinear realization. Secondly, we note that the expressions for Q\ and Q\ in 
the first line of Eq. (22) represent the leading terms in the nonrelativistic v/c expansion. It 
is not difficult to include higher order terms in this expansion as follows, 



Q\ = Ti <g) (X3 

Ql = T 2 <g> 0-3 



2\Z\ + 



2\Z\ + 



p 2 + {W'f + a 3 W" 



2M r 



p 2 + {W'Y + a 3 W" 
2M r 



1/2 



1/2 



(23) 



where the square root is to be understood as an expansion in 1/\Z\. 

In concluding this section, we note that within the context of the present Af= 2 system 
one can formulate a general statement: given the subalgebra (6) with two supercharges it is 
always possible to elevate it to a superalgebra with four supercharges and a central charge 
Z by adding the two additional supercharges (23). 



III. AN N = 2 WZ MODEL IN TWO DIMENSIONS 

A. Introducing the model 

With the aim of concretely illustrating the general arguments of the previous section, 
we now consider a specific model. A suitable example exists in two dimensions, obtained 
by dimensional reduction of a four-dimensional Wess-Zumino model with two chiral super- 
fields. The latter is the deformation of a model considered previously in Ref. [15]. The 
superpotential is 

W($,X) = — $- -$ 3 -A$X 2 + /imX 2 + — uX, (24) 
A3 A 

where $ and X are two chiral superfields, m is a mass parameter, A is the coupling constant, 
while jj, and v are deformation parameters. By an appropriate phase rotation of fields and 
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the superpotential one can always make m and A real and positive. The parameters fi and 
v are in general complex, 



fj, = /i! + ifi 2 , v = ui + iv 2 . (25) 

The four real dimensionless parameters fii, fi 2 , v \ and v 2 will form our parameter space 
{fi}. For technical reasons the parameter ji = fj, 1 + ifi 2 will be assumed to be small in what 
follows, /ii j2 <C 1. Furthermore, we will consistently work in the approximation in which the 
superpotential is linear in /i; this corresponds to terms of 0(fi 2 ) in the scalar potential. This 
limitation is not a matter of principle but, rather, for technical convenience. In the limit of 
small fi we can obtain all formulae in closed form. We will also take the coupling constant 
A to be small, A/m < 1, so that a quasiclassical treatment is applicable (except on some 
exceptional submanifolds in the parameter space). 

As a two-dimensional model, this theory has extended M= 2 supersymmetry, and ex- 
hibits solitonic kinks interpolating between the distinct vacua. In two dimensions they are 
particles (in four dimensions they would be domain walls). The dimensionality of the BPS 
supermultiplet is two, while that of the non-BPS supermultiplet is four. 

Substituting 

* = £.(U + V), X = £(U-V), (26) 
we arrive at the following action, 

S = ^ |^ J d 2 xd*6 (UU + VV) + Jd 2 xd 2 6W(U,V) + h.c)j} , (27) 
where the dimensionless superpotential W is 

1 1 

W(U, V) = U - - U 3 + V - - V s + - (U - V) 2 + v (U - V) . (28) 
3 3 2 

The vacua of the model are defined by dW/du = 0, dW /dv = 0, 

1 + v — u 2 + ji{u — v) — , 

l-u-v 2 -fi{u-v) = 0. (29) 

For real ji and v the solutions to these equations define four different vacua with real 
values of the fields and real values of the superpotential W(w, v). The vacuum structure is 
illustrated in Fig. 3 for small /x. One of these vacua (denoted as {++} in Fig. 3 ) corresponds 
to a maximum of the real function W(u, v) on the real section of the variables u and v, the 

other vacuum (denoted as { }) corresponds to a minimum of W(u, v), and the remaining 

two vacua ({H — } and { — h}) are saddle points. 

In this situation there exists [15] a continuous family of real BPS solitons, i.e. of solutions 
to the real BPS equations 

Id dW Id dW 

m dz ou m dz ov 
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interpolating between the { } vacuum with W m i n at z = — oo and the {++} vacuum 

with W max at z = oo. All these solitons are degenerate in mass: M = W max — W m i n , and 
can be viewed as a superposition of non-interacting primary solitons: one going from the 
vacuum with W m i n to one of the saddle points, and the other soliton going from the saddle 
point to the vacuum with W max - The parameter labeling the solutions in this family can 
be interpreted in terms of the distance between the basic solitons, and thus the degeneracy 
in energy implies that there is no interaction between the basic solitons at real \i and u, at 
least for some finite range of the these parameters. 



Re v 



{- 


+} 




{-1,0} 


{+ 


+} 


{0,-1} 














{0,1} 


{- 


-} 


{1,0} 


{+ -} 



Re u 



FIG. 3. Structure of vacua and solitons in the Rett, Rev plane for real v and /x . 

The decoupling of the dynamics of the primary solitons at // = is trivial, as the 
superfields U and V are also decoupled within the underlying field theory. However, at 
Hi 7^ 0, there is no such decoupling within the field theory but, nevertheless, the primary 
solitons do not interact at rest (provided [i and v are real). This is a manifestation of the 
nontrivial "no-force" condition for BPS states. 

B. Decoupled solitons, /i = case 

At fi = the model is extremely simple: the fields U and V are not coupled. Their vevs 

are 

u = ± u + , v = ± v_ , (31) 
where we introduce the notation 

u± = VTT^ . (32) 
The masses of the BPS solitons are given by 

15 



4 m 3 
3 ~X* 



n u v\ + n v u^_ 



(33) 



where the topological charges are n u>v = 0, ±1 (see Fig. 3). 

However, as noted in [15], not all combinations of charges are realized. For a generic 
value of the complex parameter v only the {1, 0} and {0, 1} solitons and their antiparticles 
exist as BPS states. To have a BPS state with both n u and n v nonvanishing, one needs to 
align in the complex plane the two terms, v\ and z/ 3 , contributing to the mass in Eq. (33). 
The relevant conditions are 



Im 




n„ 



Rc 




> . 



(34) 



These conditions define a curve in the complex v plane presented in Fig. 4. 




This curve is the curve of marginal stability for the model. In the case under considera- 
tion, with no interaction, the CMS coincides with stability domains for composite solitons, 
they only exist on this curve. 

The curve in Fig. 4 consists of three parts which can be parametrized as 

7T 

z/ = tanhcr, lmcr = 0, ±— . (35) 

3 

The part sitting on the real axis between v = ±1 (corresponding to Imcr = 0) is the 
stability domain for the {1, 1} composite solitons (and their antiparticles). The other two 
parts, Ima = ±7r/3, give the stability domain for the {1, —1} and { — 1, 1} solitons. 

The bifurcations at v — ±1 are due to the vanishing of the mass of one of the primary 
solitons at these points. It is explained by the degeneracy of vacua at these values - instead 
of four vacua only two remain at v — ±1 (strictly speaking there are still four, but they 
coalesce in pairs). These are simple analogs of the Argyres-Douglas points [13] in gauge 
theories. 
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C. Stabilization by ji 



The model at \i = is a very degenerate case. Indeed, the extra {± 1, ± 1} states exist 
only on the CMS and are nothing but systems of two noninteracting {±1,0} and{0,±l} 
solitons. The relative separation between the primary solitons is an extra classical modulus, 
on quantum level the {±1,±1} solitons are not localized states. As we will show, the 
introduction of a nonvanishing Im/x = /i 2 expands the domain of stability for the extra BPS 
state which then occupies a finite area near the original curve. Thus, setting /i 2 nonzero 
leads to an attraction of the primary solitons. 

Using n as a perturbation parameter we find the vevs and values of the superpotential 
W for the four vacua to first order in /i. 

{++}: + 

W++ = v\ + | ul + ii (1 - v-v + ) ; 

{+ - }: - = , + + f(i + ^), . = - ,- + f( 1 + ^), 
{ - +}: „=-„ ++ |( 1+ £), »=,. + f(i + £), 

>V_+ = -- i/j. + vl + n (1 + v+v-) ; 

W- = ^ - ^! + // (1 - v + v-) . (36) 

The BPS masses are given by |Wjj — WVj'| and the alignment conditions which define the 
CMS to first order in ji become (c.f. Eq. (34)), 

Im(4^r = 0, Imf4±^V /2 = , (37) 

where the conditions clearly differ only by a choice of the branch of the square root in the 
terms linear in jj,. Analytical expressions for the CMS are simpler in terms of the complex 
parameter of a (related to v by Eq. (35)). In the complex cr-plane the CMS is given by the 
curves 

<7 2 = ± yU 2 COsh — - COsh 1 / 2 <7i , 



<j 2 = — ± sinh — ^ Re 
3 2 



H cosh 1 / 2 ^<7i + i 
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^2 = -77 




) 



(38) 



where the indices 1 and 2 refer to the real and imaginary parts, o = o\ + io?,- 

The curves of marginal stability in the v plane are presented in Fig. 5. They form the 
boundaries of the stability domains for the composite BPS states marked in the figure. 



FIG. 5. The domains of stability for the composite BPS states (shown for fi2 = 0.2). The hatched 
region along the real axis is the stability domain for the {1, 1} solitons and its antiparticles, in the cross 
hatched one the {1, —1} solitons and its antiparticles are stable. 

Figure 5 exemplifies different metamorphoses of the composite BPS solitons on the CMS: 
crossing some boundaries leads to disappearance of the BPS state from the spectrum, on 
others the original BPS state disappears but a new one appears. 

The figure also shows exceptional points on the CMS, where two stability subdomains 
of the same BPS soliton touch each other. We shall address a dynamical scenario at such 
points in Sec. IV C. Note also four points of bifurcation (the Argyres-Douglas points) where 
a pair of the vacuum states collide. 



In this subsection we will find a solution to the BPS equations for the composite {1, 1} 
soliton. The construction explicitly demonstrates that in the vicinity of the CMS this soliton 
is a loosely bound state of the primary constituents. For definiteness we choose the region 
near the real i/-axis and the { } — > {++} transition. The BPS equations have the form 




D. A loosely bound composite BPS state 



1 d« 

m dz 



c 



1 + is* - (u*f + fi*(u* -v*) 



1 dv 

m dz 



c 



1-u* -(v*) 2 -n*(u* -v*) 



(39) 



where 
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We will use perturbation theory in //. The part of the CMS chosen for consideration at 
zeroth order in \i corresponds to real v. — 1 < v\ < 1, u 2 = 0. Then, at this order, a = 
and the solution for u and v reads 



u 



(0) 



ui + tanh [ui + m (z — z u )\ , — v\- tanh m (z — z v )\ , 



(41) 



where v\± = y/ui±l (see Eq. (32)) and the parameters z u and z v are arbitrary and denote 
the positions of the centers of the u- and f-solitons. 

At first order in /x, the soliton solutions become complex. With an expansion about the 
leading order solutions u^°\ v^ of the form 



u = u {0) + (ui + iu 2 ) + 
the equations (39) lead to 
1 d 



v = v (0) + ( Vl + iv 2 ) + 



(42) 



-2u^ui + m(u^ -v^) , 
Vi = —2v^ Vi — Hi (u 



m dz 
1_ d_ 

m dz 

— -7- u 2 = 2m (0) m 2 - v 2 + a 
m dz 



"<>_„((>)) ) 



h - (u^r 

2v^v 2 + u 2 + a vl_-{v^f 



1 d 

m dz 



v 2 



- A*2 
A*2 



u (0) _ v (o) 



u (0) _ v (o) 



where 



3 z/i+ - vi 
ol = - v 2 



*4+ + v \- 



(43) 



(44) 



Let us consider the equation for u 2 . The function cosh 2 [z/ 1+ m(z — z„)] is the solution of the 
homogeneous part of this equation, and the full solution is 



u 2 (z) = cosh 2 m(z — z u )\ x 



J — c 



r s —v 2 + a 
in / da; 

-00 



v 



— A*2 



-u(°)(a;) 



cosh [vi + m (x — z u )\ 



(45) 



As z — > —00 the solution satisfies the boundary condition 



lim 1x2(2) 



2z/i 



A*2 



^1- 



(46) 



consistent with Imw in Eq. (36) at the order considered here. As z — > 00 the solution 

^2(2) grows exponentially unless the relation 
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J — ( 



dx 



-v 2 + a 



uf + -(u^(x)) 



1*2 



u M(x)-vW(x) 



cosh [v 1+ m (x — z u )\ 



oo boundary condition u 2 



(47) 
lmu ++ is also 



is fulfilled. Once this relation is met the z 
satisfied. 

The relation (47) can be viewed as a constraint ensuring orthogonality of the inho- 
mogeneous part in the ^-equation and the zero mode cosh -2 [ui + m(z — z u )\ in u±. This 
approximate zero mode corresponds to a shift of the -u-soliton center and at the same time 
also represents the spatial dependence of the corresponding fermionic zero mode 2 . The re- 
lation (47) then fixes the separation z u — z v of the two primary solitons and can be presented 
in the form: 



w 



^Z u Z v ) 



1*2 K 



where 



1 

K — - 

2 



and the function w u of the soliton separation s — z u — z v is defined as 

da; 



1 f'°° dx 
Wu(s) = - / -2— tanh 

2 J-oo cosh X 



— x + z/_ m s 



cosh x 

It is important that the condition (48) also ensures that the solution for t> 2 , 
v 2 (z) = cosh 2 [ vi- m (z — z v )\ x 



(48) 



(49) 



(50) 



m 



z v 2 + a 
dx 

-oo 



(yW(x)) 2 ] +H2 \uW(x) -v^{x) 



cosh m {x — z v )\ 



(51) 



is finite at both z — > +oo and z — > — oo, and thus satisfies the proper boundary conditions. 
As for the solutions for the real parts, u± and vi, described by the first pair of equations 
in (43), these solutions always exist, due to the existence of real BPS solitons in the model 
with real parameters [15], as discussed in Sec. Ill A. Thus no additional constraint arises. 

Here we make a few remarks on the properties of the function w u (s), defined by the 
integral (50). The symmetry properties of this function can easily be seen by writing it as 
the derivative w„(s) = dg(s)/ds of the function 



9u(s) = — / dz {1 — tanh m z] tanh \y_ m (z + s)]} , 

2 J—oo 



(52) 



which is symmetric under separate reversal of the sign of v or s. Thus w v (s) is even in 
the index v. w- v (s) = w v (s), and is odd in the variable s: w u (—s) = —w„(s), and is 



2 In the next section we will show that this orthogonality condition is equivalent to the vanishing 
of a particular supercharge. 
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monotonically increasing from w u (s — > — oo) = —1 to w u (s — > +00) = +1. At large positive 
s its asymptotic behavior is given by 



w v (s) = 1 [ u + exp(— 2 v_ms) — v_ exp(— 2 v + m s)\ + 



(53) 



where the ellipses stands for higher powers and mixed products of the two exponents: 
exp(— 2 v- m s) and exp(— 2 v + ms). At v — the integral in equation (50) can be expressed 
in terms of elementary functions, 



Wq(s) = cothm s 



m s 



sinh 2 m s ' 

and the asymptotic behavior of wq(s) as s — > +00: 

w (s) = l-(4ms + 2) e~ 2ms + ... 



(54) 



(55) 



is in agreement with the v — > limit of the expression (53). Plots of w„(s) for a few values 
of v are shown in Fig. 6. 




FIG. 6. Plots of w„(s) at v — (solid), ^ = 0.8 (dashed), v = 0.95 (dot-dashed), s is measured in units 
of 1/to. 

The limited magnitude of ^(s), |^(s)| < 1, means that the BPS solution we consider 
only exists in the range 



1^2 1 < I A*2 1 « • 



(56) 



As expected the boundaries of this range coincide with the part of CMS found previously 
by algebraic means near the real v axis. 

It is then simple to find the value for the distance Sq between the primary solitons in the 
BPS composite state. Say, for v x = 0, we have when \/i 2 \ — 1^2 1 <^ | /x 2 1 



_,m I so| 



r) In 77, where 77 





A*2 




1^2 




^2| 



(57) 
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IV. QUANTUM MECHANICS OF TWO SOLITONS 



The BPS state which connects the vacua {++} and { } and exists within the stability 

domain can be viewed as a bound state of one w-soliton, located at z = z u , and one f-soliton, 
located at z — z v . The equilibrium separation between the solitons, s = z u — z v , at which 
the minimum is achieved at given u 2 and 1x2, is determined from the equation (48). In this 
section we consider the super symmetric quantum mechanics of the two soliton system. The 
SQM system describes the BPS bound state (which is the ground state in the problem) 
within the stability domain, as well as low-lying non-BPS exited states. 

The formulation of this problem refers to an effective description of the two solitons 
as heavy particles with masses M u and M v in terms of their coordinates z u and z v . This 
approximation is natural near the CMS where the binding energy is small relative to the 
soliton masses. For slowly moving solitons \z u \, \z v \ <C 1 the nonrelativistic energy can be 
written as 

E = M u i + M v Z ^ + U(s) , (58) 

where the dot denotes the time derivative, and U (s) is the interaction potential depending on 
the separation s = z u —z v between the solitons. Separating out the center of mass coordinate, 
we come to a standard quantum mechanical Hamiltonian for the relative motion, 

H-^U(s). (59) 

The supersymmetric generalization of this Hamiltonian is given in Eq. (20) which depends 
on the superpotential W(s). Below we will find the expression for this SQM superpotential 
by comparing the field theoretic supercharges evaluated on the soliton solutions with the 
SQM realization in Eq. (22). An alternative derivation of W based solely on conventional 
bosonic considerations is presented in an Appendix. 



A. SQM superpotential from field-theoretic supercharges 

The action (27) for the Wess-Zumino model leads to the following expression for the 
supercharges Q a , 

Q = ^_^ 2 J dz d t u ip + d z u 7°7 1 -?/; + im <9 S W 7*V* + (u — > v , ip — > 77) , (60) 

where u, ip a and v, r\ a are the bosonic and fermionic components of the U and V superfields, 
and W(u, v) is the superpotential of the model. The two remaining supercharges Q a are 
just the complex conjugates of Q a . 

Let us first evaluate the supercharges for the -u-soliton in the leading approximation, i.e. 
when V2 — ii\ =/i2 = 0. The field u is given by Eq. (41), u = = z/i + tanh [vi + m (z — z u )}, 
while v is a constant, v = For the fermionic fields we substitute zero modes, two of 

which are in the field ip ai and there are none in r\ a , 
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zero modes 



'ib u 
\ a u 



(61) 



In this expression M u = (4m 3 /3A 2 )z/f + is the mass of the w-soliton, a u and b u are real 
fermionic operators entering as coefficients of the normalized zero modes, and their algebra 
is fixed by canonical quantization, 



a l = h l = 1 , { a u, K} = . 
Upon these substitutions, the supercharge Q a becomes 



(62) 



(63) 



which can be rewritten in terms of real charges (see Eq. (15) in Sec. II C for definitions), 

Q\ = ^2M~ u a u , Qj = ^2M~ u b u , Q\ = Q\ = Q. (64) 

The result for the supercharges matches the general construction of Sec. II C wherein the 
operators a u and b u can be realized as Pauli matrices, e.g. a u = t\ and b u = r 2 . 

Now let us find the supercharges corresponding to the {1, 1} configuration of the u- and 
f-solitons at v 2 — ii\ =/i 2 = 0. We choose boosted soliton solutions, 



u 



(0) 



= v 1+ tanh 



V\_ tanh 



u 1+ m\z — z, 



V 



v\_ m z 



V 



(65) 



where p is their relative momentum, and the total momentum is zero. The fermions are 
given by Eq. (61) for ip and by a similar expression for i] with the substitution u — > v, where 
u- and -y-fermions anticommute. With time- dependent solutions the terms d t uip, d t vr) now 
contribute to the supercharges (60), 



Q{ = J2(M U + M V ) (a u cos5 + a v smd) , Q\ = J2(M U + M V ) (b u cos 5 + b v sin 5) , 



p 



{—a u sin 5 + a v cos 5) , 



Ql 



p 



{—b u sin 5 + b v cos 5) , 



(66) 



where we have defined cos 5 = \J M u / (M u + M v ) . We observe that the relative motion 
implies that the 'composite' state is non-BPS in the absence of any interaction between the 
solitons. 

In order to switch on the interaction we consider nonzero fi and To obtain the 
result to first order in these parameters it is enough to substitute the same leading order 
expressions for the bosonic and fermionic fields accounting for the terms linear in fi and v 2 
in the expression (60) for the supercharges, as well as for the phase a of the central charge. 
The linear dependence on \x and v 2 arises from the terms 
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m 



V2X 2 



dz 



- %v 2 7 (ip* - rf) 



(67) 



in (60). The phase a is also linear in v 2 (see Eq. (44)), and needs to be taken into account 
in Eq. (15) when relating Q a with Q]f . 

The resulting supercharges are (Q\, Q\ are not changed and are written here for com- 
pleteness) 



Q\ = p{M u + M v )a, 
1 



where we denote 



pa + W\s) b 



Ql = y / 2(M M + M v )b, 
1 



Qi 



pb-W'(s) a 



(68) 



a = a u cos 5 + a v sin 5 , 
a = —a u sin 5 + a v cos 8 



b = b u cos 5 + b v sin 5 , 
6 = — b u sin 5 + b v cos 5 . 



(69) 



The quantum-mechanical superpotential (more precisely its derivative W) is then read 
from Ql (or Q\) to be 



W\s) 



3M r 



[ll 2 KW ul (s) -V 2 \ , 



where 



M r . 



M U M V 
M u + M v 



2 m 3 (1 

3 ¥~ 



K 3 + + <4 3 -) 



(70) 



(71) 



and the function w u (s) is defined by Eq. (50). 
The SQM Hamiltonian then has the form 



HsQM = (Q 2 ) — (Ql) 



2M r 



p 2 + {W\s)f -iW"(s)ab 



(72) 



An explicit matrix realization of the four operators a u>v and b U)V , satisfying the Clifford 
algebra, can a priori be chosen in the factorized form: a u — T\ <g) 03, b u = r 2 <g> a 3 , a v — I<S>(Ji, 
and b v = I <S> cr 2 . This factorized form of the fermionic operators realizes a description in 
terms of two independent particles. This choice is perfectly acceptable and realizes the Af— 
2 superalgebra (18). However, it differs from the specific realization (22) by an orthogonal 
rotation of angle 5. In order to match the conventions used in Eq. (22) for a description 
of the two-soliton system, one has to use an equivalent representation of these operators, 
obtained by the inverse rotation: 



a u = Ti<S> cr 3 cos 5 — I <g> U\ sin 5 , b u = r 2 <g> <r 3 cos 5 — I <E> <y 2 sin S , 
a v — Ti <S> C3 sin 5 + / <S> <J\ cos 5 , 6^ = r 2 <8> a 3 sin <5 + / <8> a 2 cos 5 . 



(73) 



The final expression for the full quantum Hamiltonian of the two-soliton system can be 
written as 
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p 2 9M r [li 2 kw u1 (s) - v 2 f 3 li 2 kw' u1 (s) 
HsQM= 2M~ r + -2 (T^!y + 2 {l-ul) * 3 (74) 

with M r given in Eq. (71) and we use the matrix representation (73) for the fermions 
(omitting the tensor product with unity in H S q M ). 

It is worth noting a couple of limits in which the potential simplifies, and can be expressed 
in terms of elementary functions. Recall first of all that when v\ = the function Wq(s) is 
known analytically (see Eq. (54)). If v\ is not too large, i.e. v\ < 0.5, there also exists a 
convenient simplified form in which the superpotential is very closely approximated by the 
expression 

^a Ppr ox(s) = 3M r [/i 2 tanh(m S ) - v 2 \ , (75) 

where the reduced mass M r is taken at v\ = 0. In this superpotential we recognize the 
simplified model discussed in Sec. II B (see Eq. (7)). Another simple case arises when v\ is 
close to 1, i.e. 1 - i/J < 1, 

W'= 3Mr (fi 2 ms-u 2 ). (76) 

The potential in this case reduces to that of the harmonic oscillator. 

In the limit of large separation, the potential energy in the Hamiltonian (74) tends to a 
constant which depends in the sign of s, 

tt uf jl ^ 9Mr (±A*2«- v 2 ? 

£4 = tf(*-±oo) = - {x _ v2)2 , (77) 

(Note, however, that the spin dependent <7 3 term does not contribute). These constants 
denote the energy levels at which the continuum states appear while the ground state, 
which is the {1, 1} BPS soliton, is a zero energy eigenfunction of i^sQM- 

The origin of the two continuum thresholds is that at nonzero fx the classification for 
solitons we introduced at ji = is no longer sufficient — the w-soliton interpolating between 
the { — h} and {++} vacua (see Fig. 3) is different from the M-soliton interpolating between 

the { } and {H — } vacua, and a similar distinction arises between the v- and -5-solitons. 

Thus, the system under consideration at large s describes two channels: the u plus v solitons 
at positive s, and the u plus v at negative s. It is straightforward to verify this by calculating 
the two binding energies, 

TT? 

AE + = M M — M u — M v — — [ - - - - - W__ | ] , 

(78) 

771 

AE_ = M M — + M e = — [ - - - >V__| - - W+_| ] , 

A 

from which we observe that AE± = —U±. Note that, although the quantities AE± are of 
second order in /i 2 and v 2 , it is sufficient to use the expressions (36) which are only valid 
to first order in ji (and are formally exact in v) for the values of W^. This is due to the 
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fact that for real v and [i the values of Wy are real and AE± vanishes. Thus AE± arises as 
an effect quadratic in the imaginary parts of the differences of which by themselves are 
linear in v 2 and /i 2 . 

Finally, we also write down the asymptotic behavior of the potential as s — > ± oo, 

U(s) — ► U± + K±exp(-2v 1 _m\s\) + ---, (79) 
where the coefficient of the leading exponential term is 
3 (u 1+ + I/i_)^ 2 « 



■6M r (/i 2 KT^)+mff3(l-^>i- ■ (80) 



We have made the assumption here that v\_ < We see that the characteristic distance 
s is defined by l/mvi- which as expected is the wavelength for the lightest particle in the 
model. We also observe that the spin dependent term contributes to the exponential tail. 
Moreover, on the CMS where /i 2 k; =F ^2 = 0, it is the only contribution. This 'fermionic 
dominance' takes place in a very narrow region near the CMS, 

(1 — v?) Vi_ m 

I^T-frK 6 K- (81) 

The effect of this regime of enhanced quantum corrections near the CMS will be considered 
in more detail in the next subsection. 



B. Properties of the two-soliton system 

As expected, the second term in the potential in Eq. (74) is of higher order than the 
first in the loop expansion parameter A. However the second term is of lower order in the 
small parameter /x 2 and, by tuning /i 2 , one can study this potential both in the classical 
limit corresponding to /i 2 A 2 /m 2 and in the quantum limit /i 2 «C A 2 /m 2 , or anywhere in 
between as long as the condition for validity of the formula (74), /i 2 1, is maintained. 
Upon a slightly more detailed inspection of classical vs quantum effects in the Hamiltonian 
(74) one can readily see that in fact the quasiclassical parameter in this system is not just 
the ratio A 2 / (m 2 /i 2 ) but is determined by the parameter 

_ W»»(s ) (g2) 

introduced in Sec. II, where so is the classical equilibrium separation determined by (48). 
Recall that £ measures the quantum correction to the curvature of the potential near the 
classical minimum, the system being essentially classical for £ <C 1, and highly quantum for 

£»l- 

For the model at hand we find 



<(ao) A 2 



2m/2 2 yjl - v\ (w' Vi (s Q )y 



1 ■ (83) 
m 2 
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Near the CMS the equilibrium distance s is large, and £ takes the form 



K 



A 2 



(84) 



where for defmiteness we have again assumed that v\_ < Notice that the condition 

|£| 3> 1 agrees with Eq. (81) which, as discussed above, defines the essentially quantum 
regime in the narrow region along the CMS. 

When the system admits a supersymmetric ground state, the corresponding wave func- 
tion ipo(s) can always be found as 



ipo(s) = const exp [— W(s)] = const exp 



Mr H2K9v i. (s) ~ v 2 ms 
m 



(85) 



where for defmiteness we again assume /i 2 > 0, and g u (s), defined by Eq. (52), is the integral 
of w v (s). Independently of the quasiclassical parameter £ the maximum of ipo(s) is always 
located at s . However the spread of the wave function, i.e. the dispersion of the distance 
between the solitons in the BPS bound state, essentially depends on the parameter £. As 
£ — > the full potential has a minimum at s — sq, and the system is classically located at 
the minimum. At larger £ the minimum of the full potential shifts towards s = 0, reaching 
s = in the limit £ ^> 1, but the maximum of the wave function is still at s = sq. In the 
latter extreme quantum limit the system resembles the deuteron: the wave function spreads 
over distances much larger than the size of the interaction region. In the two-soliton system 
this behavior is even more drastic at large £ than in the deuteron: the wave function reaches 
its maximum far beyond the interaction region. The classical and the quantum behavior of 
the system at different values of £ is illustrated by a series of plots in Fig. 7, 8. 
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FIG. 7. Plots of the full potential U(s) (arbitrary units) at v\ = 0, ^2/^2=0.95 for several values of £. 
The classical equilibrium point is at ms «2.56 and is shown by heavy dot. Fig. a: details of the potential 
near minimum for £ = 0.009 (solid), £ = 0.9 (dashed), and for £ = 4.45 (dot-dashed). Fig. b: the potential 
shown at a larger scale. The curves for £ < 4.45 are practically unresolvable and coincide with the solid 
curve, the dashed curve corresponds to £ = 125. It can be noticed that the latter value of £ still corresponds 
to moderate values of A/to: A 2 /m 2 w 7.7^2- 
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FIG. 8. Plots of the ground state wave function tpo(s) for £ = 17.7 (dashed) and £ = 177 (solid). As 
above, these parameters correspond to ^2/^2 = 0.95 and the classical equilibrium point is at ms s=a 2.56. 

One may also note that in general the interaction of the two solitons is strong only at 
distances of order m~ l near s = 0: the potential is asymmetric in s and changes rapidly near 
s = 0, i.e. the force is strongest when the solitons substantially overlap. In a narrow region 
near the CMS, given by Eq. (81), an attraction at short distances creates an essentially 
quantum state, resembling a deuteron. Deeper into the stability region an exponentially 
shallow minimum of the potential at large Sq results in a quasiclassical bound state. 

Once one crosses the CMS, the wave function (85) is no longer normalizable, and the 
physical ground state of the system is non-supersymmetric. The broadening of the wave 
function for the bound state near the CMS is exhibited in Fig. 8. Thus, as discussed in 
Sec. II B the bound state level reaches the continuum on the CMS, where it completely 
delocalizes, and on crossing the CMS the {1, 1} bound state is no longer present in the 
physical spectrum. 



C. Another dynamical regime: Extra Moduli on the CMS 

The scenario discussed above, involving a short range superpotential which remains finite 
on the CMS, is only a generic description for the near CMS dynamics in certain systems. As 
discussed in Section 2, a different dynamical scenario is possible if there exist extra moduli 
on the CMS. However, it turns out that the two-field model also exhibits a dynamical regime 
of this type, and we observe in this case that the approach to the CMS is still characterized 
by derealization of the bound state wavefunction, albeit in a somewhat different manner to 
the case considered above. 

Firstly, recall that in the example considered above with [i 2 7^ 0, the approach to the 
CMS was determined by Eq. (48). In the 'interior' domain, |z/ 2 | < |/i2 1 k, the equation 
W'(s) = has a solution and consequently the composite soliton was BPS saturated. Upon 
approach to the CMS, the zero of W'(s) runs to infinity and, after crossing the CMS at 
1 1/2 1 > |//2|«, there is no longer a solution to W'(s) = and hence no BPS soliton. This 
scenario is illustrated by see Fig. 9a.. 

Now we consider a different dynamical regime, see Fig. 9b, where, in both the 'interior' 
and 'exterior' regions, the spectrum of BPS states is the same (although possibly rearranged). 
In this case one still has derealization on the CMS one still has derealization, although 
only for the wave function in this case as there is no diverging (classical) separation of the 
constituents. 
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FIG. 9. Possible scenarios for the BPS spectrum, taken from a small region of Fig. 5, where Pi and 
P2 refer to the two primary solitons, while C refers to the composite kink: (a) A composite bound state C 
exists only on one side of the CMS; (b) a bound state exits on both sides of the CMS. 

To this end let us set v = (i.e. discard the term linear in U, V in the superpotential 
(24)). As explained in Sec. Ill A, at v — the CMS is very simple, 

Im/i = ji 2 = 0. (86) 

The SQM system (A. 10) one arrives at in this case is described by the superpotential, 

TfL 

W > (s)=2 — f , 2 w (s), (87) 

where wo(s) was defined in (54). For our purposes it is important that u>o(0) = 0, and that 
there are no other zeros of wo(s). Thus the solitons always overlap classically. However, as 
one approaches the CMS the wave function still spreads out due to the fact that /i 2 — > 0. In 
particular, at large \s\, 

w {s) — ► sgn(s), (88) 

and one observes that the zero energy bound state exists for both positive and negative /z 2 
(see Fig. 9b). The wave function at large s is 

2-^-//2S0"3 1 (89) 

times either | J,) or | f) depending on the eigenvalue of 03. As \/i2\ — > the bound state level 
approaches the continuum spectrum while the wave function swells. At /x 2 = the wave 
function is completely delocalized and there is no binding. 

As alluded to above, this dynamical regime is distinct from that considered previously 
where W remained finite on the CMS; rather the CMS was characterized by the escape of 
the root of W'(s) = to infinity. In contrast, in the example considered here the root of 
the equation W'(s) = does not shift at all. Despite this one may note that £ still diverges 
near the CMS due to its inverse dependence on /i 2 - 

In fact, precisely on the CMS the potential vanishes, and thus a new quantum modulus 
arises corresponding to the relative separation of the constituents. 



29 



V. DYONS IN SU(3) M = 2 SYM 



We turn now to consider similar phenomena in M— 2 SYM. To study a model which 
exhibits a CMS in the weak coupling region, one approach is to extend the gauge group to 
rank greater than one 3 . Here we consider one of the simplest examples of this kind with 
gauge group SU(3). In the Coulomb phase this theory exhibits BPS dyon solutions with 
electric and magnetic charges associated with either of the unbroken U(l)'s of the Cartan 
torus. After choosing a convenient basis of simple roots for the algebra, one can classify the 
BPS monopole solutions into one of two types: those whose magnetic charge is aligned along 
a simple root - 'fundamental monopoles' - and those whose magnetic charge is aligned along 
the non-simple root. These 'composite monopoles' generically possess CMS curves at weak 
coupling, and so their dynamics in this regime is amenable to a semi-classical consideration. 

Composite dyons in this, and the closely related Af = 4 system, have recently been 
studied in some detail [16,18-22,25], with the conclusion that the low energy dynamics 
of two fundamental dyons at generic points of the Coulomb branch acquires an additional 
potential term. This term is associated with the misalignment of the adjoint Higgs vevs of 
the two dyons, and leads to the formation of composite BPS dyons as bound states in this 
system. We will review some of these results below, and emphasize the implications for the 
dynamics in the near CMS region. The removal of the composite state on the CMS again 
arises through derealization. 

A. The BPS mass formula 

We first review the features of the classical BPS mass formula for J\f— 2 SYM with higher 
rank gauge groups (see e.g. [22,26]), limiting ourselves to SU(N). For the consideration of 
solitonic mass bounds, we need consider only the bosonic Hamiltonian which has the form, 

H = 2Trj d 3 x {^{Ei) 2 + ^(B t ) 2 + A)$* D $ + A$* D& + ^ [<£>*, $] 2 | , (90) 

where and Bi (i = 1,2,3) are the electric and magnetic fields, and $ is the complex 
adjoint scalar ($ = ($x + i& 2 )/V2 in terms of the two real adjoint scalars). We use the 
normalization Tr T a T b = (l/2)5 ah for the generators. 

The classical vacua satisfy [$*, $] = 0, thus requiring the vev of $ to lie in the Cartan 
subalgebra H, 

($) = 0-H. (91) 

Note that the remaining Weyl freedom may be fixed by demanding that Re <fi • fl a > for a 
given set of simple roots {/3 a }. This defines a region which for |</>-/3 a | 3> A coincides with the 



Alternatively, one can add hypermultiplet matter with a large mass. In this case there is a 
discontinuity in the spectrum of quark monopole bound states on a CMS curve, which has been 
studied by Henningson [24], and the mechanism involves derealization in a manner analogous to 
that discussed in this section. 



30 



semiclassical moduli space of the theory. In this region we can safely neglect field-theoretic 
perturbative and nonperturbative quantum effects. We will only consider the case where 
the gauge group is maximally broken to U(l) Ar ~ 1 , for gauge group SU(iV). 
For a soliton solution we may define the charge vector Q 

Q ■ H = (q + ig) ■ H = f dS l (E t + iB{) , (92) 

where use of the unitary gauge is implied. The real (q) and imaginary (g) parts of each 
component of Q have the interpretation of electric and magnetic charges in the corresponding 
U(l); they are are quantized and form a lattice spanned by simple roots, 

Air 

q = q a f3 a = en*(3 a , g = g a (3 a = — nf (3 a . (93) 

e 

Here e is the gauge coupling, and and are the integral electric and magnetic quantum 
numbers. We also normalize the simple roots f3 a with the conventions (/3 a ) 2 = 1, /3 a±1 • (3 a = 
— 1/2, so that the coroots coincide with the roots. 

For general bosonic configurations, there is a Bogomol'nyi mass bound following from 
(90) which takes the form, 

M>M&x\Z±\, Z+ = V2<t>*-Q, Z_ = V2cj)-Q, (94) 

This bound is saturated by solutions of the Bogomol'nyi equations, 

Bi = Dib, E i = D t a, (95) 

along with the equation Dfa =F e 2 [6, [b, a\] = which, making use of (95), expresses Gauss' 
law in the gauge Aq = =)=a. The fields a, b are real and imaginary parts of exp(ia) $ where 
the angle of rotation a is defined in terms of the charges Q (see e.g. [16,22,27]). 

In the framework of the J\f= A supersymmetry algebra the parameters Z± are realized 
as central charges, and it is advantageous to view the system in this context (implying six 
instead of two real scalars $). Within J\f—A SUSY it is generally the case that ^ |-2L|, 
and states which saturate the Bogomol'nyi bound M > Max \Z±\ will preserve only four 
of the sixteen supercharges, and will thus be 1/4 supersymmetric. If, however, \Z + \ = 
\Z_\, states which saturate this bound will preserve 1/2 of the A/"=4 supersymmetry. The 
possibility of having 1/4 BPS states, which only occurs for gauge groups of rank larger than 
one, dramatically increases the number of CMS curves accessible to semiclassical analysis, 
since 1/4 BPS states generically exhibit regions in the parameter space where they become 
marginally stable with respect to "decay" into 1/2 BPS states. In this sense it is useful to 
think of 1/4 BPS configurations as composite. 

The discussion above was framed within A/ = 4 SYM, but this was simply for orientation. 
In order to preserve any fraction of supersymmetry, four of the six real adjoint scalars must 
vanish asymptotically, and thus the configurations discussed above all "descend" to give 
classical solutions in Af— 2 SYM. The difference is that now only Z_ remains as a central 
charge [22] and all states saturating the bound M = \Z_\ are 1/2 BPS states from the 
point of view of the M— 2 SUSY algebra. As noted in [22] , those charge sectors for which 
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\Z-\ < \Z + \ will have no BPS states from the point of view of the A/ = 2 system. Indeed, in 
this case states with M = \Z_\ are not allowed because M > \Z + \. Thus \Z_\ > \Z + \ is a 
necessary condition for the existence of Af=2 BPS states. 

Restricting our attention now to the gauge group SU(3), then on the Coulomb branch 
the gauge group is broken down to the Cartan subalgebra U(l) 2 and there can be field con- 
figurations which are electrically and magnetically charged under either of these U(l)'s [28]. 
Following Weinberg [28] we use the term 'fundamental dyons' to refer to those configurations 
whose charges are aligned along one of these simple roots. Configurations whose charges are 
aligned along non-simple roots (i.e. /3 1 + f3 2 ) will be referred to as 'composite'. We shall 
focus on a particular composite configuration which has received considerable attention in 
the recent literature - namely the composite dyon which has equal magnetic, n a M = (1, 1), 
and differing electric, n% = (<?i/e, g2/e), charges along the simple roots. 

B. Marginal stability and Coulomb-like interaction 

The BPS mass formula for the (1,1) dyon with n a M = (1, 1), n% = (qi/e,q2/e) takes the 
form, 

M (1)1) = \Z_\ = V2\( qi + ig) (f> ■ (3 l + (q 2 + ig) <f> ■ f3 2 \ , (96) 

where g = 4-7r/e. This configuration has a CMS curve where the (1,1) dyon is marginally 
stable with respect to two fundamental dyons: the first is aligned along (3 1 , and has n a M = 
(1,0), n% = (qi/e,0), while the second is aligned along /3 2 , and has n a M = (0,1), n% = 
(0, (72/e). The masses of the fundamental dyons are 

M a = V2 I (q a + tg ) ( f>.f3 a \=m a (l + ^ + ) , (97) 

where in the second equality we have made use of the fact that e 2 <^ 1 in order to write the 
dyon mass as a small correction to the mass of the corresponding fundamental monopole, 
m a = V2g\4> ■ (3 a \. Introducing 00 a as the argument of the vevs, 

4)-(3 a = \(f>- f3 a \e iu)a (98) 

we see that the marginal stability condition M( ljl ) = M 1 + M 2 fixes the argument of the 
ratio 

UJ — U)\ — Cl>2 , (99) 

to be equal to the argument of the ratio of complex charges Q2/Q1 where Qi = qi + ig. This 
implies that the CMS equation is u = uj c where oj c is defined as 

sinc c = , {q2 - q ] )9 . (100) 
vV + <7i 2 \A 2 + 5 

Provided n a E and n a M are of a similar order, the angle oo c is small in the limit e 2 < 1, i.e. 
when the electric corrections to the dyon masses are much smaller than the corresponding 
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monopole mass as in (97). Thus in this limit the vevs <fi ■ (3 a are only slightly disaligned, 
and we can make use of an expanded version of Eq. (100), 

co c = ^ (l + 0(q*/g 2 )) , Aq = q 2 - qi . (101) 

We are now in a position to verify the general claim of Sec. II A that the Coulombic 
interaction vanishes on the CMS. At large distances dyons can be viewed as point charges 
which interact at rest through electrostatic, magnetostatic, and scalar exchange. The elec- 
trostatic and magnetostatic interactions are fixed by the corresponding charges, while the 
scalar exchange can be read off from the asymptotic form of the Higgs field of one of the 
primary dyons (in a physical gauge where the configuration is a linear superposition of the 
fundamental dyon solutions [27,29]). The effective Coulombic interaction then takes the 
form, 



Coul 



8nr 



(102) 



Similar expressions have appeared in [27] and [19]. One observes that on the CMS, where 
the angle uj is given by Eq. (100), the Coulombic potential Vc ou i vanishes. 
When expanded to second order in qi/g, the potential Vc ou i takes the form 

where Aq = q 2 — q\ is defined in Eq. (101). If the vevs were aligned, i.e. uj = 0, we see that 
to quadratic order, the potential is repulsive [29,30] (as opposed to the SU(2) case [31]) and 
depends only on the electric charge difference Aq. However, for (guj/Aq) > 1 the potential is 
attractive and the (1,1) BPS dyon exists with a mass M^ i) given in the same approximation 
by 

M (1>1) = \Z_\ « (Mi + M 2 ) - M r (Ag ~r )2 > > i) > ( 104 ) 

where M r = mim 2 / '(mi +m 2 ) is the reduced mass of the monopoles, and the corresponding 
masses Mi )2 and mi^ are defined in Eq. (97). On the other side of the CMS in the range 
\gu/Aq\ < 1 we have repulsion and the the (1,1) BPS dyon does not exist. 

Its interesting to note that the Coulombic potential reverts to attractive form once more 
when (gu/Aq) < — 1, and the (1,1) bound states reappear. This has a simple interpretation 
in the framework of Af—4 supersymmetry: the lowest mass in this range saturates the \Z + \ 
central charge (note that now \Z + \ > \Z^\), 

M (M) = \Z + \ « (Mi + M 2 ) - M r (Ag + r )2 , < -l) . (105) 

In terms of A/ = 4 SUSY the (1,1) state at (gu/Aq) < —1 is 1/4 supersymmetric, but 
preserves a different subalgebra as compared to the (gu/Aq) > 1 case above. Moreover, the 
generators of this subalgebra are not part of the M= 2 superalgebra. In terms of M= 2 this 
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means that the supermultiplet is not shortened, but nonetheless the Bogomol'nyi bound is 
saturated at the classical level 4 . Thus, we see an interesting example where the "BPS" 
nature of the state does not imply multiplet shortening. The presence of these states in 
M=2 SYM was also noted in [22]. 

We now address the question of what happens to the (1,1) state on the CMS, i.e. when 
\gu/Aq\ = 1 and the 1/r terms in the potential vanish? As we will see in Sec. VD the dy- 
namics on the CMS is governed by repulsive 1/r 2 terms demonstrating, even at the classical 
level, that there is no localized bound state on the CMS. 



C. Zero Modes and Moduli Spaces 

We will shortly consider the low energy dynamics of the fundamental dyons comprising 
the (1,1) system. However, we first recall a few details regarding the zero mode structure 
of dyon solutions in M= 2 SYM. For BPS dyons in pure M= 2 SYM the unbroken M= 
1 supersymmetry is enough in this case to pair the bosonic and fermionic zero modes [32] 
so we shall focus here just on the bosonic modes. Generic dyon solutions, corresponding 
to the embedding of the SU(2) monopole along some root of SU(3) have four bosonic zero 
modes [26,28] parametrizing the moduli space, 

Mi = R 3 x S 1 . (106) 

These modes are naturally identified as the center of mass position in R 3 and the S* 1 is an 
isometry conjugate to the conserved electric charge. 

For dyons embedded along a simple root, this is the moduli space for all choices of 
field-theoretic moduli. However, if we consider composite monopoles, then the monopole 
moduli space M. enlarges to a space of dimension eight, as is compatible with separating 
the constituents into two isolated fundamental monopoles [26,27]. This result was obtained 
in [26] using the index calculations of Weinberg [28] for real Higgs fields. 

For the case at hand, the magnetic charge is g = g{f3 l + f3 2 ) and asymptotically the 
eight dimensional moduli space is simply M. 1 x M. x . However, its exact form has also been 
deduced in [29,33-35], 

^, = R »x^22t. (107) 

The first R 3 factor corresponds to the center of mass position, while the second R factor 
refers to the coordinate conjugate to the total electric charge. The corresponding metric is 
flat. The relative moduli space Mtn is positive mass Taub-NUT space (which is asymptot- 
ically R 3 x S* 1 ). Its four coordinates describe the relative distance r between the cores, 
with the corresponding polar and azimuthal angles 9 and </?, and also the relative phase \i 
conjugate to the relative electric charge Aq. The factor Z denotes a discrete identification 



^his saturation will be lifted by field-theoretic quantum corrections. 
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for the charge coordinates, ensuring that the asymptotic geometry has a compact factor 
S* 1 x S 1 , associated with the conserved charges. 

The Taub-NUT metric, in our conventions 5 takes the form, 



d<4 N = g ^ dz"dz u = m(r) dr 2 + r 2 d9 2 + r 2 sin 2 6 dip 



e 2 m(r) 



dx~\ — cos 9 dip 
2 



with a 'running' mass parameter, 



m(r) = M r + — , 



e 2 r 



, (108) 



(109) 



which asymptotes to the reduced mass M r when the relative separation r diverges. 

In terms of the internal U(l) angles of the fundamental monopoles, the combination 
£ = ifti + tp2 is conjugate to the total electric charge q t (or more precisely, to qt/e), 



m x q\ + m 2 q2 
Qt = ; , 



(110) 



while x — (mi^2 — ^2^1)/ ( m i + m 2) is conjugate to Aq/e, i.e. to the relative electric charge 
introduced above [34,35]. 



D. Moduli space dynamics 

As first discussed in this context by Manton [36] , the low-energy dynamics of fundamental 
monopoles may be understood as geodesic motion on the underlying moduli space. This 
picture extends to dyon solutions with aligned charges, but recent work on the dynamics of 
the fundamental constituents of the (1,1) system [16-22] has shown that for two fundamental 
dyons with misaligned charges the Lagrangian following from the geodesic approximation 
needs to be corrected by a new term [19,22]. In this subsection we will partially review these 
results, while emphasizing the features of the near CMS region. 

The construction of Refs. [19,22] can be reformulated in terms of the following Lagrangian 
for the relative moduli z^ = {r, x} 

L I * = \g™* l ? + g™* t c?', (in) 

where the metric g™ is given by Eq. (108) and the 'potential' G v , which is a Killing vector 
generating the x isometry, is 

G v = -M r u5". (112) 

9 x 



5 We follow [22] with the exception that % is rescaled to have a period of 27r rather than Att, and 
consequently the conjugate momentum is integer ((92— 5i)/e) rather than half- integer ((52 — <Zi)/2e) 
valued. 
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Here oo is the angle of disalignment between the condensates 4> ■ (3 a , see Eq. (98). The term 
containing C 1 is dynamically significant due to nontrivial fibering of the S l associated with 
X in the Taub-NUT metric. 

The classical Hamiltonian then has the form 

#rci = \ 9tn K - Gjfa - G„) , (113) 

where = g™{z u + G") are the canonical momenta. In terms of the original field theory 
this Hamiltonian is interpreted as M— \Z_\, and thus BPS states are 'vacuum' states of this 
Hamiltonian. 

Two of the canonical momenta, namely n v and n x , are conserved quantities conjugate to 
isometries along the azimuthal angle tp and the phase x- The value of is the ^-projection 
of the angular momentum l z , and n x is equal to Aq/e. Substituting these and the inverse 
metric (j^N into Eq. (113) we obtain 

H «* = o~tt + + o ( ?2- 2 fl - t- cosd ) 2 + v{r) ' (114) 

2m(r) 2m(r)r 2 2 m(r) r 2 sin 9 \ 2e / 

where 

v ^-^( 1+ ^)\ Aq -'^ +Aq ^h)' (115) 

is the only term in H IC \ which depends on the field-theoretic moduli (via a;). 

This Hamiltonian vanishes at -n r = 0, n e = 0, and the equilibrium values r and 9 of the 
corresponding coordinates are given by 

2-k Aq e 

r o = ~YT1 a ' cos9 = 2l z — , (116) 

e z M r goo — Aq Aq 

when (guo/Aq) > 1 and \2l z e/Aq\ < 1. There is no solution for (guo/Aq) < 1, i.e. the 
BPS state ceases to exist upon crossing the CMS where gou — Aq. We see that the sys- 
tem describes the composite state as a bound state of the dynamics whose spatial size, 
corresponding to the separation of the primary constituents, diverges on approach to the 
CMS. 

It is instructive to expand the potential V(r) at large r 

Wr)=i^ <*i^ + <Mztogg + (£gg + (117) 

where we omitted 1/r 3 terms and higher powers of 1/r. The constant term V(r — > oo) 
marks the start of the continuum. Indeed, adding \Z_\ from Eq. (104) we obtain Mi + M 2 
in the limit r — > oo. The 1/r term coincides with the Coulombic potential (103) discussed 
earlier. It provides attraction for \gu/Aq\ > 1, the range where the bound states exist. The 
range {guo/Aq) < —1 corresponds, as discussed above, to Mnu = \Z + \. What we see in 
addition is the repulsive 1/r 2 term which leads to the existence of an equilibrium at large 
r near the CMS. However, on the CMS it becomes the dominant term, and so there is no 
localized state on the CMS. 
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We conclude this section with some brief comments about quantization. The quantization 
of the dyon system in the context of Af— 2 SYM was first discussed in detail by Gauntlett [32] 
in the case where the Higgs vevs are aligned, and this discussion has since been generalized 
to the case considered here by Bak et al. [20] for Af= 4 , and by Gauntlett et al. [22] for 
Af— 2 SYM. The crucial feature of this system is that the relative moduli space inherits a 
triplet of complex structures j( a \ a = 1..3, and is a hyperkahler manifold. Consequently, 
the system exhibits Af— 4 supersymmetric quantum mechanics with four real supercharges 
Qa where the index A can be associated loosely with a quaternionic structure = (1, J a ). 
These supercharges satisfy the superalgebra, 

{Qa,Qb} = 25 A bH iA , (118) 

where H re \ is the supersymmetric completion of the Hamiltonian defined in (113), to be 
interpreted as M — \Z-\. One can compare this with the SQM constructed in Sec. 2. 

An interesting feature of this system is that the symmetries of the superalgebra and 
the moduli space combine to ensure that the wavefunctions have a nontrivial dependence 
on the angular moduli, as well as the relative separation r. Specifically, the ground state 
wavefunction has the functional form [20] \& = ^o(r,a a ) where a a are the basis 1-forms 
on the S 3 parametrized by (9, tp, x)- This dependence is hinted at through the ^-dependent 
terms in the Hamiltonian (114) above. One may speculate that because of the high degree 
of symmetry in this system - the bosonic system possesses an additional conserved quantity 
of Runge-Lenz type [25] - a more precise separation of variables may be possible, but we will 
not pursue this issue here. We note only that, as demonstrated above, the derealization on 
the CMS is associated with the cancellation of the terms of 0(1/ r) and depends purely on 
the relative separation. Moreover, this conclusion in the classical bosonic system apparently 
extends to SQM [20]. 

VI. DELOCALIZATION VIA MASSLESS FIELDS 

On particular submanifolds of the CMS, the discussion we have presented above may be 
incomplete because certain fields may become massless. Indeed, generically there will be 
particular points on the CMS where states which are stable on both sides are massless. The 
presence of these singularities in moduli space can then be thought of as the 'origin' of the 
CMS since marginally stable states may not be single valued on traversing a contour around 
the singularity, and so a discontinuity in the spectrum becomes necessary for consistency. 
For example, this point of view provided one of the first arguments for the fact that the W 
boson must be removed from the spectrum inside the strong coupling CMS in Af— 2 SYM [3] . 

In this section we will discuss the behavior of BPS states near these singular points. 
Within a simple 2D Wess-Zumino model we will find that the discontinuity of the BPS 
spectrum is explained by the derealization of fermionic zero modes of the soliton on the 
CMS. The CMS in this case corresponds to the 'collision' of two vacua in the parameter 
space, and thus one might anticipate similar phenomena in Af — 2 SYM near Argyres- 
Doulgas points [13] when the singularities associated with monopole and quark vacua collide. 
Unfortunately, this occurs at strong coupling and is out of the range of our semi-classical 
analysis. 
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A. Breaking N = 1 to N = 1 and the restructuring of WZ solitons 



The model is a simplified version of the M= 2 Wess-Zumino model [14] in two dimen- 
sions considered in Sec. III. We shall set the second field X to zero, but consider a new 
perturbation which breaks M= 2 down to A/ = 1 . This setup was introduced in Ref . [37] (see 
Sec. 8). 

The superpotential prior to perturbation is then of the standard Landau- Ginzburg form, 
and its worthwhile recalling a few pertinent details of these theories. A general classification 
of the Af = 2 Landau-Ginzburg-type theories in two dimensions was given in [1], while 
construction of the representations of the J\f = 2 superalgebra with central charges was 
presented in Refs. [23]. It was shown that the supermultiplet of BPS soliton states is 
shortened, and this shortened multiplet consists of two states {u, d} as we discussed earlier 
in Sec. III. In particular, in J\f—2 theories there exists a conserved fermion charge /. The 
fermion charge of the u and d states is fractional but the difference is unity, f u — fa — 1. 

What changes on passing to A/ = 1 in two dimensions? The irreducible representation of 
the M= 1 algebra for BPS states is now one-dimensional (to the best of our knowledge this 
was first noted in Ref. [37]). The only remnant of the fermion charge is a discrete subgroup 
Z 2 which is spontaneously broken. 

It is natural then to expect a restructuring of the BPS spectrum when M= 2 is broken 
down to J\f= 1. We will study the manner in which restructuring occurs by considering 
the spectrum of fermionic zero modes of a soliton solution as we vary the soft breaking 
parameter /z. We observe that for small /x the BPS spectrum remains the same as in the 
unbroken Af=2 theory. However, starting at a critical value /z* - corresponding to a 'point of 
marginal stability' - half the BPS states disappear from the spectrum. This occurs because 
quasiclassically the counting of states in the supermultiplet is related to the counting of zero 
modes of the soliton and when fi reaches /i* some of the fermionic zero modes become non- 
normalizable. To follow their fate one can introduce a large box. Then the number of states 
does not change, but at /i = the identification of states with zero modes implies that half 
the BPS states spread out all over the box while for they lie on the boundary of 

the box and are removed from the physical Hilbert space. This picture is quite analogous 
to the quantum mechanical discussion in Sec. II. However, as we shall see, the quantum 
description is complicated here by the presence of a massless field. 

We take the Kahler metric to be canonical and the cubic superpotential of the model is 
conveniently represented in terms of real bosonic variables <pi, % = 1, 2, 

th? A 

W(y?i,y9 2 ) = — ^--y^ + Ay^ . (119) 
In fact, W(y?i,y?2) is harmonic 

d 2 W 



dtpidtpi 



= for AT = 2. (120) 



as it is the imaginary part of the four dimensional superpotential which is analytic in ipi+iifi2- 
This is therefore a reflection of M= 2 supersymmetry in two dimensions. 

Now, to break A/=2 down to M= 1 we consider a more general, nonharmonic, superpo- 
tential W (ipi, lfi2), 
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W(y?i,y2 2 ) = — (ft - -(pl + Xtp^l + ^tpl, (121) 
where jj, is the soft breaking parameter. There are two vacuum branches, 



r V" 



<* = ±2V <* = °j' = ^ = ±J ^T- ' (122) 



but the second exists only for fi > m, and vacua collide when /j, — m. 

This model exhibits a classical kink solution which interpolates between the first set of 
vacua, and is given by 

Tfl TYIZ 

<Pi = 2\ tanh "2"' ^ = °- ^ 123 ^ 



(124) 



It satisfies the classical BPS equations, 

d Vi _ dW 

dz dipi 

The zero modes corresponding to this kink are as follows: 

(a) one bosonic mode: 

Xo = C 1 (125) 
cosh (mz/2) 

in the field ifi corresponds to (the spontaneous breaking of) translational invariance, 

Xo oc d^i/dz. 

The constant C in Eq. (125) is a normalization constant; its explicit numerical value is not 
important. (Below the normalization constants in the zero modes will be omitted.) 

(b) two fermionic modes: The first zero mode in the field ipi^ (the indices number the su- 
perfields and fermionic components in the basis where 7 = 02, 7 1 = icr^) has the same form 
Xo a s the translational mode. It is not accidental, the corresponding differential operators 
are the same due to A/= 1 supersymmetry. The second fermionic zero mode 

_ exp(-//z) 

6 " cosh 2 (m^/2) (126) 

appears in the field ip 2 ,i- At /i = the existence of this mode is a consequence of the A/ = 2 
SUSY, and £0 coincides with xo- At nonvanishing //, when the extended SUSY is broken, 
this zero mode is maintained by virtue of the Jackiw-Rebbi index theorem [38]. 

An interesting feature of the zero mode £ is that it is asymmetric in z for /j ^ 0. 
Moreover, this mode is normalizable only for 

fi<m. (127) 

This is readily seen from its asymptotics, 
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£o(z — » oo) ~ e 



-(p+m)z 



oo) 



~(n-m)z 



(128) 



The explicit form of the zero mode for few values of \i is exemplified in Fig. 10. 



-10 




5 
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FIG. 10. The kink profile (solid line), with the zero mode £o for /V m = (dotted line), /z/m = 0.8 
(short-dashed line), and fi/m = 1 (long-dashed line). Note that the vertical scale has been altered for ease 
of presentation. 

The loss of normalizability occurs at n — m, when 



in one of the vacua between which the soliton solution interpolates. In other words, one of 
the vacuum states has gapless excitations at this point. Indeed, in the z — > — oo vacuum, 
the eigenvalues of the fermion mass matrix are: m and fj, — m, and thus indeed at the point 
(1 = 171 where the vacuum branches meet, there is a massless field. This system represents 
a simplified analog of an Argyres-Douglas point in that the massless field arises through 
the collision of vacua. Furthermore, we see that this infrared effect destabilizes one of the 
fermionic zero modes of the soliton. 

To study the infrared behavior in detail let us put the system in a large box, i.e. impose 
boundary conditions at z — ±L/2 where L is large but finite. We choose these conditions 
in a form which preserves the remnant supersymmetry in the soliton background (i.e. the 
BPS nature of the soliton), 



(see Ref. [37] for details). It is easy to check that the soliton solution (123) as well as the 
zero modes (125), (126) are not deformed by these boundary conditions, i.e. (123) remains 
a solution of the classical BPS equations with the appropriate boundary conditions at finite 
L. 






(129) 
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In the finite box there is no problem with normalization; the zero mode (126) remains a 
solution of the Dirac equation in the soliton background for all //. However, at // > m the 
mode is localized on the left wall of the box instead of sitting on the soliton as is the case at 
ji < m. Thus, at // = m the critical phenomenon of derealization starts. As we will show 
below, upon quantization this means that some BPS soliton states have disappeared from 
the physical Hilbert space. 

B. Quantization 

We shall not present a detailed analysis of the quantization of the system here as it 
requires a somewhat different treatment to the supersymmetric quantum mechanics we have 
considered thus far. In this case, one needs to consider the dynamics of the light field in 
addition to the collective coordinates of the soliton. 

However, provided we only consider a region somewhat away from the CMS, the spec- 
trum is easily determined. As usual, the remnant M— 1 supersymmetry pairs the nonzero 
modes (one bosonic to two fermionic) around the soliton (see e.g. Sec. 3G of [37]), and the 
relevant contributions cancel. Thus the soliton spectrum is determined by the zero modes, 
corresponding to which we have one bosonic collective coordinate zq corresponding to the 
center of the kink, and two real Grassmann collective coordinates a\ and a.2 determined 
by the zero modes, = aiXo(z) + ' ' ' an d ^2,1 = a 2^o( z ) + ' ' '• Combining them into 
one complex parameter 77 = (ai + ia^/v^, the collective coordinate dynamics at ji — is 
determined by the quantum mechanical system 

C eff = -M + \mz\ + iMffq, (130) 

where M is the physical kink mass, which we can set to one. 

The quantization is carried out in the standard manner. If z and 77 are the canonical 
coordinates, we introduce the canonical momenta 

7^0 = ^0, TT v = -ifj, (131) 

and impose the (anti) commutation relations 

[t*o» z ] = -i , {tc v , 77} = -% . (132) 

One then proceeds to construct the raising and lowering operators in the standard manner. 
From Eqs. (131) and (132) it is clear that r\ can be viewed as the lowering operator, while 
fj is the raising operator. One then defines the 'vacuum state' in the kink sector by the 
condition that it is annihilated by 77, 

77! 'vac') = . 

The application of fj produces a state which is degenerate with the vacuum state. |'vac') 
and 7/| 'vac') are two quantum states which form a (shortened) representation of M = 
2 supersymmetry (at fj, — 0). It is clear that these two states, which are degenerate in 
mass, have fermion numbers differing by unity. 
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What happens at /i 7^ 0? At // < m the situation is exactly the same as at fi = 
(apart from the fact that the fermion number is not conserved now and we must classify the 
states with respect to their Z 2 properties). We have two degenerate quantum states, both 
are spatially localized and belong to the physical sector of the Hilbert space. At /i — m 
only the vacuum state is localized. The spatial structure of the second state fj\ 'vac') has 
a flat component, which extends to the boundaries of the box. At fi > m this component 
is peaked at the boundary. The easiest way to see this is to introduce an external source 
coupled to Eijipiipj. Thus, the state 77] 'vac') disappears from the physical sector of the Hilbert 
space. The supercharge Qi acting on the state |'vac') produces this state itself, rather than 
another state. (We recall that Q 2 annihilates |'vac').) Formally this looks like a spontaneous 
breaking of the remnant supersymmetry. 

VII. CONCLUDING REMARKS 

Using quasiclassical methods we have argued that the underlying dynamics of margin- 
ally stable solitons is determined (generically) by non-relativistic supersymmetric quantum 
mechanics. Composite BPS states which disappear on the CMS were found to do so through 
a process of derealization in coordinate space. Within the quantum mechanical description 
this process was associated with the bound state level reaching the continuum, while further 
progress beyond the CMS leads to a potential with a non-supersymmetric ground state. 
This is a generic picture. In certain cases the CMS can be a boundary between sectors with 
different composite solitons, the quantum mechanical potential then vanishes on the CMS. 

One of the crucial features allowing a detailed investigation of the effective quantum 
mechanical dynamics in the two field model considered in Sections III and IV was the linear 
realization of supersymmetry in the two-soliton sector of the non-relativistic SQM system. 
This embedding in 1+1D is similar to the 3+1 D effective dynamics of two BPS dyons in 
Af=2 and A/"=4 supersymmetric gauge theories [19-22]. In 3+1D the effective Coulombic 
interaction governs dynamics near the CMS. We demonstrated this in Section V for the 
composite dyon in SU(3), showing how the BPS states swell upon the approach to the CMS. 

One may wonder whether some of the conclusions noted above regarding the behavior 
of composite bound states near the CMS might not be artifacts of the quasiclassical ap- 
proximation. In particular, returning to pure M— 2 SYM with gauge group SU(2), the 
classic example of marginal stability with which we started this discussion is that of the W 
bosons on a CMS curve at strong coupling [4,5], for which our methods are not directly 
applicable. We are going to dwell on this issue in a separate publication [39]. Here we will 
briefly present two suggestive arguments pointing to the conclusion that this phenomenon 
involves derealization in the same manner as the examples we have discussed. 

The first observation involves duality. If we consider a region very close to the CMS for 
the W boson and not too far from the monopole singularity, we can consider a point particle 
approximation for the monopoles and dyons within the dual magnetic description. Provided 
we are close enough to the CMS, a nonrelativisitic approximation will be reliable, and from 
this viewpoint the dissociation of the W is superficially quite similar to that of the bound 
states of dyons discussed in Sec. V, with the roles of electric and magnetic charge reversed. 

The second observation involves the realization of the BPS states considered here in 
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terms of string junctions [6, 7, 16, 40] in Type IIB string theory and its extension to F- 
theory [41]. Although there are still subtleties with this realization, specifically concerning a 
mismatch between field- and string-theoretic counting of bosonic moduli [17], its interesting 
that the disappearance of marginally stable states in this framework appears to universally 
implies derealization. The crucial point is that this process involves shrinking one or more 
of the spokes of the junction to zero length, while it has recently been pointed out [42] 
that the equilibrium separation of the two constituent states in the field theory is inversely 
proportional to the length of the shrinking prongs. 

As a final remark, its worth commenting on additional subtleties which arise when consid- 
ering extended BPS objects. In particular, although we concentrated here on BPS particles, 
the notion of marginal stability is more general as supersymmetry algebras may also ad- 
mit central charges supported by extended BPS objects such as strings and domain walls. 
Indeed, our classical analysis in Sec. Ill may be lifted to four dimensions where the kink so- 
lutions describe BPS domain walls. However, we concentrated on particle states specifically 
for the reason that quantization in this case leads to quantum mechanics, which is of course 
well-understood. The main technical difficulty in extending these arguments to solitons such 
as domain walls is that in addition to the dynamics of relative collective coordinates, one 
also needs to consider the massless sector of the field theory on the worldvolume of the 
soliton. 
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APPENDIX: CLASSICAL SOLITON POTENTIAL 

In this appendix we present an alternative derivation of the classical potential (W'(s)) 2 
entering the SQM Hamiltonian. This approach is purely bosonic, but requires knowledge of 
the composite soliton solutions obtained in Sec. III. 

We start from the expression (79) for the binding energy AE + , 

AE + = -— 12 ^ 2 > . A.l 

(For z/ 2 //i 2 < the relative sign between z/ 2 and /i 2 in this expression must be reversed.) 
The formula (A.l) gives the minimum of U(s) — U+, where U + is the value of potential at 
s = oo, while the position of the minimum in s is given by Eq. (48). We can combine these 
two results in order to find the classical expression for U (s) by using the standard Legendre 
transform approach. We introduce a source term in the original superpotential (24), thus 
replacing W(<1>, X) by 
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m 



X; j) = X) -yjT, (A.2) 

where j is a dimensionless (and in general complex) parameter corresponding to the strength 
of the source. For a static classical configuration described by this superpotential the calcu- 
lation of the energy in fact gives the minimum of the functional €(j): 

£(j)=E(j)+^±Jdzd 2 6X(x,6)+h.c^j , (A.3) 

where E(j) is the value of the original energy on the configuration which extremizes the 
action at a given source strength j, and X(x,9) is the X superfield evaluated on that 
configuration. (In fact, being static, X does not depend on time.) 

Clearly the effect of the source term is equivalent to a shift in v: v — > v + j and for 
our purposes it is sufficient to consider a purely imaginary source, j = ij2- Then the s 
dependent part of the functional £{j) for the two-soliton static configuration is read directly 
from equation (A.l) after replacing v 2 by v 2 + j: 

and the relation between the j dependent equilibrium position s and the value of j is derived 
from Eq. (48) 

^2 + h = A*2 « w Vl (m s) . (A. 5) 

The quantity of interest for us here, however, is not the functional S(j) as a function of 
j, but rather the binding energy AE as a function of s. The latter is found in the standard 
way from the relation 

AE = £(j)-j^ (A.6) 

with the variable j being eliminated in favor of s, using the relation (A. 5). Performing this 
simple operation on the expressions (A. 4) and (A. 5) one finds 

AE(s) = U(s) -U + = ^4 3( " 2 ~ /Wl(mS))2 - K 3( " 2 ~ /i2/t)2 , (A.7) 

A 2 KV 1+ Vi_ X 2 KV 1+ Vi_ 

which represents the classical interaction energy of two primary solitons separated by a 
distance s. Naturally, the minimum of AE(s), as found from this expression, coincides with 
that given by equation (A.l) at the separation s determined by (48). 

Comparing the last term on the right hand side of Eq. (A.7) with the expression (77) for 
U + we see that they coincide. Thus, the potential is 

A 2 KV-l+V-l- 
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where the subscript reminds us that this is the potential found at the classical level. This 
result correctly reproduces the energy difference between the asymptotic states at both 
infinities. With this normalization, one may readily check that the the minimum (zero) of 
the potential corresponds in (79) to the mass of the BPS {1,1} soliton, 

Comparing the above expression for the potential with the general form of the SQM 
Hamiltonian, 

+ Kwf+TS''* (A ' 10) 

we readily derive the superpotential (up to a sign) 

This coincides with the result (70) derived in Sec. IV A, by evaluating the field-theoretic 
supercharges. 
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